Realization of the 7r-state in junctions formed by multi-band superconductors with a 

spin-density-wave. 
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Using a simple model for Fe-based pnictides, we calculate the Josephson current Ij in a tunnel junction 
composed by such superconductors. We employ the tunnel Hamiltonian method and quasiclassical Green's 
functions. We study both the case of coexistence of the superconducting (A) and magnetic (SDW — spin den- 
sity wave) order parameters and the case when only the superconducting order parameter exists. We analyze 
ideal and nonideal nesting and show that the current /j depends on the mutual orientation of magnetization 
of the SDW only in the case of nonideal nesting. 

It is found that the realization of the -junction is possible in both cases. We show also that in S++ /I/S+- 
tunnel junctions an anomalous term ~ cos may arise in the Josephson current and this result remains 
valid in the absence of the SDW (S++ and S+- are superconductors with (A, A) respectively (A, -A) order 
parameters in different bands and I denotes an insulating layer). 

PACS numbers: 74.45. +c, 74.50.+r, 75.70.Cn, 74.20.Rp 



Introduction. The superconducting order parameter (OP) 
in the BCS model A is either constant or angle-dependent 
as it occurs in anisotropic superconductors^. Super- 
conductors of new types discovered in last decades are 
characterized by a more complicated OP. For example, it 
is proposed that the OP in superconductors with heavy 
fermions may have the OP corresponding to the triplet 
pairing and having a vector structure^. The OP in 
high-Tc superconductors depends on angles and in some 
of them corresponds to the singlet, so-called d-vfave pair- 
ing: A(a) - Ao(cos(fcj;ajc) - cosikyOy)), where kx - fcpcoso:, 
ky - fcp sin a are the components of the Fermi momen- 
tum fcp in the CuO layers^. This means that the OP turns 
to zero in certain directions and therefore nodes arise in the 
excitation spectra. 

The OP in the recently discovered superconduct- 
ing materials — Fe-based pnictides, which have rather 
high Tc (~ 60K)'^ — is also nontrivial. These materials be- 
long to a class of multi-band superconductors; their band 
structure consists of electron and hole bands. The OP in 
these bands may have not only different amplitudes but 
also opposite signs. If the OP in the hole and electron bands 
have different signs, one speaks of the s+_ -pairing — in 
contrast to the the s++ -pairing in case of the same signs of 
the OP in different bands (see, for example, reviews 5-11). 
Investigating the structure of the OP in different types of 
superconductors is a very important task since this study 
may shed light on the mechanism of superconductivity in 
these materials. 

One of the effective methods to determine the structure 
of the OP is the measurement of the Josephson current. 
For example, the sign change of the OP in high- Tc super- 
conductors has been proven in experiments in which the 
critical Josephson current Ic was measured in a setup con- 
taining two Josephson junctions connected by a supercon- 
ducting loop (i.e., on SQUID)^. The dc Josephson effect 
in multi-band superconductors has been studied theoreti- 



cally in many papers'^^"^^. The calculations in Refs. 13 and 
14 are focused on the multi-band superconductor MgBg, 
whereas the main attention of the authors of Refs. 15-19 
is paid to the Fe-based pnictides. Agterberg et al.^^ have 
shown that the Josephson S/I/Smb junctions may have a 
negative critical current 7c if the OP A is negative in some 
of the bands (here, S and Smb mean single-band and multi- 
band superconductors, respectively, I stands for an insulat- 
ing layer) . This idea allows a simple physical interpretation. 
The Josephson current 7c in an S/I/Smb junction can be writ- 
ten as 7c oc Y.a ^^alRa, where Ra is the resistance for elec- 
tron transitions from the S superconductor to the band a in 
the Smb superconductor. It is clear that if Aq is negative, for 
instance in the band with a-1, the current 7c may be also 
negative. This happens provided the resistance 7?i is suffi- 
ciently small. The ground state of the Josephson junction 
with negative 7c is called the :^-state. 

Note that the existence of the :rt-state of the Josephson 
junction is interesting in itself because such junctions can 
be used in practical applications (see, e.g., Ref. 20 and ref- 
erences therein). The jr-state is realized in S/F/S Joseph- 
son junctions and is being studied very actively (for reviews 
see Refs. 21-25). Therefore, there is a need both from the 
point of view of fundamental research and of the future ap- 
plications to study the possible realizations of the jr-state 
in Fe-based pnictides. In all publications mentioned above 
the presence of the spin density wave (SDW) in these ma- 
terials is ignored. On the other hand, it is known that there 
is a region on the T-x-plane (temperature and doping level) 
where the superconducting and magnetic (SDW) phases co- 
exist. 

In the present paper we calculate the Josephson cur- 
rent 7c in tunnel junctions formed by multi-band supercon- 
ductors, i.e., in the Smb/I/Smb junctions, both with an SDW 
and in the absence of the SDW. We show that the critical 
current 7c in junctions with the SDW consists of two terms. 
The first one is proportional to AiAr and does not depend 
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on the angle between the magnetization vectors mj ^ in the 
SDW on the left (1) and on the right (r). The second term 
is proportional to AiAr(mi • rrir). This means that this com- 
ponent can be negative and, as we will show, it can prevail 
so that 7c oc cos(20), where 26 is the angle between the vec- 
tors mi and nir. Thus, 7c may be negative in such junctions. 
We show also, that an anomalous term I^ncosq) arises in 
S++/I/S+- junctions in addition to the ordinary Josephson 
current 7an sin (p. 

System under consideration. Model. We consider a tun- 
nel Smb/I/Smb junction. Each superconductor on the left 
and on the right is described by the Hamiltonian j which 
contains the superconducting and magnetic energies taken 
in the mean field approximation^^"^^. Transitions of elec- 
trons between superconductors is described by the tunnel- 
ing Hamiltonian 

E i-^«)3fla,raAl + H-C.} = E{C'^HTCl + H.C.}, (1) 
p.a.p P 

where the matrix elements describe the electron tun- 
neling between the same bands ifa-f}, 3~aa = STa, and be- 
tween different bands if a 7^ /3. In the latter case one has 
STaf,- ST^^ for the identical superconductors at the left and 
right. We assume that the matrix elements ,5^^ do not de- 
pend on momentum p. The band a- \ (resp. 2) is assumed 
to be the hole (resp. electron) band. The C = Cans opera- 
tors are related to the cia operators as follows: Cans - ^ns 
for a = 1 and Cans - Bns for a = 2. In the hole band one has 
Ain] - ^Jj for n = 1 and Ain\ - for n-2.ln the electron 
band the relations B2n\ - for n- \ and Bzn] - d^2\ ^^'^ 
n-2 take place. One can see that the labels a, n and s are 
the band, Gor'kov-Nambu and spin indices, respectively. 

The operator Hx can be written in terms of the matrices 
p, f and a which operate in the band, Gor'kov-Nambu and 
spin spaces: Hx = T ■ f 3. Here, the matrix f is given by 

t ^^[3-^%Qo+ sr.-i[r x2io-r' x22o]\ , (2) 

where B'+ = [5\ ± ^2) 12 and Y = y' + iV" = B\2. The ma- 
trices Xans are defined as Xans- Pa-fn-^s- It is worth 
making an important note. As is known, in tunnel junc- 
tions composed by single band superconductors, the rela- 
tion 9'{p,p') - !T* {-p,-p'] holds which is a consequence 
of the time reversal symmetry. Therefore, the matrix ele- 
ments that describe tunneling between identical bands 
are real if we assume that these matrix elements do not 
depend on momenta p and p' . However, the matrix ele- 
ment y = Sr\2 describes tunneling between different bands. 
In this case, the initial and final states are different and the 
idea about time reversal symmetry is not applicable. There- 
fore, the matrix element T , generally speaking, is complex. 

One can obtain the Eilenberger-like equation for quasi- 
classical Green's functions in the, e.g., right electrode with 
the self-energy part £x,r which is related to the tunneling 
Hamiltonian: ix.r^ f -gi-f^^". The matrix gi in the self- 
energy part £x,r is the quasiclassical Green's function in the 
left electrode. It can be found from the Eilenberger equa- 
tion neglecting the self-energy £x,r since we assume small 



tunneling probability (the method of quasiclassical Green's 
functions is described in reviews 33-37). In the case of ideal 
nesting, this matrix function is given in Ref. 30. If the nest- 
ing is not ideal, this function acquires a more complicated 
form. 

Josephson current. From the generalized Eilenberger 
equation one obtains the rate of the charge variation with 
time, e.g., in the right electrode: (<3fQr -1- (3(/Qr)U=f' = 7x, 
where Qr - eN^{Q) J deTr{X3oog^} with the density of states 
at the Fermi level N-^ (0) and gf is the Keldysh component of 
the Green's function. In equilibrium, this function is equal 
to gf = (gj^-g/^)tanh(ej3) with /3 = {2T)~^. The tunneling 
current in equilibrium is the nondissipative Josephson cur- 
rent 7j. It is given by 

7j = ci(4OTr)^Tr{X33o[fgi(w)f,gr(w)]}, (3) 

(J) 

where ci = 7reA/i(0)A/r(0)/16 and gi,r(w) are the Green's func- 
tions in the left (right) electrodes in the Matsubara repre- 
sentation. 

Ideal Nesting. Using Eqs. (2-3) we calculate the Joseph- 
son current between superconductors like Fe-based pnic- 
tides separated by a tunnel barrier. We start with the sim- 
plest case when the SDW is absent. Then, the trace in Eq. (3) 
is not zero only for the condensate component which has 
the form (we assume lAhoiel - lAeiectronl - 

^ A [l323Cos((^/2) + Xoi3sin((p/2) for 5+_ -pairing, 

® Xo23 cos((p/2) + X313 sin((/?/2) for s++ -pairing, 

(4) 

where Q> - vZj^ + A^ and q) is the phase difference between 
the left and the right superconductors. We calculate the cur- 
rent 7j for different junctions: a) S++/I/S++, b) S+-/I/S+- 
and c) S++/I/S+-. In the symmetric cases of identical su- 
perconductors forming the junction we obtain the standard 
formula 7j - 7c sincp with different critical current 

2jl-^il' + l-^2l' + 2«(r2) forS++/I/S++, 
Jc/7o-MI 1 15^-^12 + |2-25R(y2) forS+_/I/S+_, ^' 

where 7o = [nM2eRn) tanh(A/2r), = {S'lf + {S'lj + \r\^ 
and R-^ ^ 47ie^NimNA0) {\^\l + \^\l + \y\^] is the resis- 
tance of the junction in the normal state. In the case of an 
asymmetrical S++/I/S+_ junction we obtain a quite differ- 
ent result 

7j/7o = |^|-2[(|^i|2-|^2l^]sing> + !3(y2)cos(^] . (6) 

As follows from Eqs. (5) and (6), in case of real tunnel- 
ing matrix elements ^12 and T the critical current may be 
negative in junctions S+-/I/S+- and S++/I/S+-. In the first 
case, 7c is negative if the interband transitions dominate, 
i.e., =3^ 2 « y- In the S++/I/S+- junction, 7c changes sign 
at < \32f. These results resemble those obtained in 
Refs. 13-19, where the Josephson current in junctions of the 
type S/I/Smb has been studied. 

However, if the tunneling matrix elements are complex, 
the obtained results do not reduce to those established ear- 
lier^^"^^. Especially interesting is the result for the critical 
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FIG. 1. (Color online.) Considered setup of the Josephson junction. 



current in an S++/I/S+- junction. In this case, the Joseph- 
son current has the form: Ic oc 7ci sm{2(p) + Id cos(2(^). This 
means that a spontaneous current arises even at zero phase 
difference; or a finite phase difference is established across 
the junction in a disconnected circuit. The presence of a 
spontaneous condensate current has been established ear- 
lier in different systems where the time-reversal symme- 
try breaking takes place. For example, this current arises 
in a superconducting loop containing a :7r-Josephson junc- 
^jQj^38,39 gp biiayer*". The presence of the sponta- 

neous current can be understood from the physical point 
of view as follows. Even if the global phase difference is zero 
[2q) - 0), the phase difference between the band with nega- 
tive A and the band with positive A is not zero. Due to this, 
the interband transitions lead to the nonzero component 
current 7j. Note also, that there is a similarity of this effect 
and the anomalous proximity effect in S+- system studied 
in Refs. 41 and 42. 

Nonideal nesting. Now, we consider the most interest- 
ing case of the non-ideal nesting when the superconduct- 
ing (A) and magnetic (SDW) order parameter may coex- 
ist in a certain interval of doping level x and tempera- 
ture T^^^^, i.e., we introduce a parameter Sfip, describing 
the mismatch of the effective Fermi surfaces of the bands. 
Then, the Green's function for the 5+_ -pairing, in case cp-0 
and the magnetization vector oriented along the z-axis, has 
the form g+_ = g+_(0,0) = goso^oso + gioo-^ioo + gusXus + 
g2i3^2i3 + g300^300 + g323^323- For the 5++ -pairing we have 

g++ = g++ (0, 0) = g023^023 + g030^030 + gl23-f 123 + gl30-f 130 + 

g2i3-^2i3 + g300-^300- That is, the matrices g+_ and g++ have 
a rather complicated form. In our notations for the S+- sys- 
tem we follow previous papers^^'^" and the quantities with 
tilde denote the according ones in the S++ system. However, 
only two terms— g323l323, gioo-^ioo in g+- and go23^023, 
gi30^i30 in g++ — contribute to the Josephson current. 

If the phases (pi 2 are different {+q)/2] and the angle 6^0 
(see Fig. 1), the Green's functions in the left (right) su- 
perconductors g\jiq),6) are expressed through the matri- 
ces g(0,0) with the help of the unitary transformations: 
g\,r{(p,d) = R±e ■ S±^ ■ g{0,0) ■ Rig ■ §1^, where the signs + re- 
late to the left (right) electrodes. The transformation ma- 
trices are: S+^p - exp (+1X330(^/4) and R+q - exp(+iX33i0/2). 
They can be called the rotation matrices in the Gor'kov- 
Nambu and spin spaces. 



In the case of S++/I/S++ and S+-/I/S+- junctions we ob- 
tain for the Josephson current 

i^/iooism(p-\3-\-^Y^{g'ol3[m^+m^+2my^)] (?) 

+g[locosmm^2-\yf)} 

and 

Ij/Io<x sm(p-\^\-^Y.{S323[\^i\^ + \^2\^ -2my^)] (8) 

+g[locosmm^2-\yf)} , 

respectively where g^23 = ^(+^\V\~^'^{{^mo ~ 
Sfil + g'uQ = w„AM/MoCVr^9W. g'wo = 

MVMoC'ixr^^ix] and g'^^ ^ ^ Ac'ixr^mx] + 

5/Xp!3{x}), with yf = ^Wlj^ + A2 + a)l-Sfil-2(+, 



and ( - \J A^ + o)^. We see that the angle-dependent part in 
Eqs. 7 and 8 differs from zero only if Sfip ^ (the condition 
of the non-ideal nesting). 

In the case of the S++ /I/S+- junctions we obtain no angle 
dependent part 

7j/7ocx \3-\-^Y.{s'oi3g'o23[\^i\^ -\^2\^]s^^V + ^iy^)coscp} . 

(9) 

This formula looks like the corresponding expression for the 
case of an ideal nesting. 

71 -state. In principle, the realization of the :rt-state is possi- 
ble in all three cases. The case c) of the S++ /I/S+- junction is 
discussed above. In the cases a) andb), i.e., in the S++/I/S++ 
and S+_ /I/S+_ junctions it is possible to have the change of 
the sign of the critical current near Tc for sufficient large ra- 
tio of Ts/Tc. 

However, the :^-state is most easily realized and observed 
in the S+-/I/S+- junction — at arbitrary temperature. To 
achieve that, one needs to adjust the tunneling probabili- 
ties related to the matrix elements S'ap in a way as to make 
the first term in Eq. 8 (proportional to gj^j) vanish. As- 
suming the matrix elements 3'a real, this can be achieved, 
e.g., for the following ratios: ^2/^1 - 1, 5R{^i2}/^i - V2, 
9{5l2}/!ft{^i2} = In this case, die critical current 

7c oc cos(20) is proportional to cosine of the angle be- 
tween the magnetization vectors in the leads. This, to- 
gether with the mechanically favorable properties of the Fe- 
pnictides^^'^^, makes the S+-/I/S+- junction very attractive 
for possible applications in the quantum devices. In partic- 
ular, the possibility of creating wires allows one to think of 
the realization of the so-called ^-junction (arbitrary phase 
shift, not only or 71) out of pieces of Fe-pnictide wire, put 
together mutually rotated and separated by an insulating 
layer — a chain of Josephson junctions. 

Discussion. Starting from the generally accepted model 
for the Fe-based pnictides and a generic tunneling Hamil- 
tonian, we investigated the Josephson current in the junc- 
tions consisting of two Fe-pnictides separated by an insu- 
lating layer: a) S++/I/S++, b) S+-/I/S+- and c) S++/I/S+-, 



4 



where the indices indicate the pairing of the superconduct- 
ing OP of the system. We find that, even in case of perfect 
nesting, the critical Josephson current shows unusual prop- 
erty for the S++/I/S+- junction being finite for zero global 
phase difference (the so-called 0-junction). In case of per- 
fect nesting, there also exists the possibility of a sign-change 
in the S+-/I/S+- and S++/I/S+- systems, thus enabling a 
:/r- state in the junction. 

However, most easily the n^-state is realized in the case of 
non-ideal nesting in the S+- /I/S+- junction, where it is pos- 
sible to achieve a cosine dependence of the critical Joseph- 
son current on the mutual orientation of the magnetiza- 



tion vectors in the superconductors, 7c cx cos (20). Impor- 
tant ingredient for this effect is the co-existence of the SC 
and the SDW phases, as it occurs in the Fe-pnictides. It is 
not possible to achieve the pure cosine dependence of Ic on 
the mutual orientation of the magnetization vectors in the 
S++/I/S++ junction because the pre-factor of g^^g in Eq. 7 
is always finite. In the S++/I/S+- junction, there is no de- 
pendence on the mutual orientation of the magnetization 
vectors of the SDW. 

The authors are grateful to I. Eremin for useful remarks 
and discussions. 
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